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G^ ; ABSTRACT 



We analyze the dynamics of 2-D stationary, line-driven winds from accretion disks in cat- 
f^ ' aclysmic variable stars. The driving force is that of line radiation pressure, in the formalism 

'^ ■ developed by Castor, Abbott, & Klein for O stars. Our main assumption is that wind helical 

^vj , streamlines lie on straight cones. 

(^ ' We find that the Euler equation for the disk wind has two eigenvalues, the mass loss rate and 

^\ ■ the flow tilt angle with the disk. Both are calculated self-consistently. The wind is characterized 

by two distinct regions, an outer wind launched beyond four white dwarf radii from the rotation 
axis, and an mnerwind launched within this radius. The inner wind is very steep, up to 80° with 
the disk plane, while the outer wind has a typical tilt of 60°. In both cases the ray dispersion is 
O I small. We, therefore, confirm the bi-conical geometry of disk winds as suggested by observations 

*-^ ' and kinematical modeling. The wind coUimation angle appears to be robust and depends only on 

Ci • the disk temperature stratification. The flow critical points lie high above the disk for the inner 

wind, but close to the disk photosphere for the outer wind. Comparison with existing kinematical 
and dynamical models is provided. Mass loss rates from the disk as well as wind velocity laws 
are discussed in a subsequent paper. 

Subject headings: accretion disks — cataclysmic variables — hydrodynamics — stars: mass-loss 
— stars: winds 



1. Introduction 

Accretion disks are ubiquitous in astrophysical systems ranging from new-born stars to compact objects, 
like white dwarfs, neutron stars and black holes, both stellar and galactic. Due to their high temperatures 
and large surface areas, disks appear to be among the most luminous objects in the Universe. Strong 
dissipative processes which accompany accretion around compact objects can eleviate the radiation energy 
density in and above the disk, leading naturally to radiation-driven winds, similar to winds from hot stars. 
Observational signatures of such winds have been unambigously detected in cataclysmic variables (CVs) 
(Heap et al. 1978; Krautter et al. 1981; Klare et al. 1982; Cordova & Mason 1982) and in active galactic 
nuclei, hereafter AGNs (Arav, Shlosman, & Weymann 1997, and refs. therein), but their understanding 



proved to be challenging for theorists. In this and the following paper (Feldnieier, Shlosman, & Vitello 1999; 
hereafter Paper II) we focus on different aspects of disk winds in CVs, such as their 2-D geometry, solution 
topology, mass loss rates and velocity profiles. AGN disks will be discussed elsewhere. 

Theoretical understanding of winds from accretion disks is hampered by their intrinsically multi- 
dimensional character and by the richness of various physical processes supplementing the basic hydro- 
dynamics of the flow. A number of different driving mechanisms for disk winds have been predicted and 
analyzed, from magnetic torques to X-ray disk irradiation (i.e., Compton- heated and thermally-driven winds) 
to resonance line pressure (e.g., Blandford & Payne 1982; Begelman, McKee, & Shields 1983; Cordova & 
Mason 1985; Woods et al. 1996). Disks in non-magnetic CVs with high accretion rates, ^ 10^^ M© yr^^, 
have an energy output which peaks in the (far-)ultraviolet, similarly to O, B and WR stars. Their spectra 
exhibit features which bear similarity to those found in hot and massive stars, and which are attributed 
to winds driven by radiation pressure in resonance and subordinate lines of abundant chemical elements, 
so-called line-driven winds (LDWs). Observational evidence in favor of LDWs winds from hot stars and 
disks includes but is not limited to the P Cygni line profiles of C iv, N v and Si iv, ionization levels, high 
terminal velocities and their correlation with the luminosity, and UV line behavior during continuum eclipse 
in CVs. 

Pioneering works by Lucy & Solomon (1970), Castor (1974) and Castor, Abbott, & Klein (1975; CAK 
in the following) have shown that O star winds result from scattering of radiation in the resonance lines 
of abundant elements. Elegantly formulated theory of the LDWs from O stars by CAK, Cassinelli (1979), 
Abbott (1980, 1982), Pauldrach, Puis, & Kudritzki (1986) and others (for a textbook account, see Lamers & 
Cassinelli 1999), was successfully applied to individual objects. Further refinements of this theory by Owocki 
& Rybicki (1984, 1985) and Owocki, Castor, & Rybicki (1988) addressed the issue of stability of the flow. 

First application of the LDWs to accretion disks (Shlosman, Vitello, & Shaviv 1985; Vitello & Shlosman 
1988) emphasized the non-spherical ionizing continuum and driving force as well as a bi-conical geometry 
of the outflow. Under a broad range of conditions disk atmospheres in CVs and AGNs become dynamically 
unstable because the line opacity effectively brings them into a super-Eddington regime. Continuum photons 
absorbed by the UV resonance lines and re-emitted isotropically contribute to the momentum transfer to 
the wind. This process can be described as a resonant scattering which conserves the number of photons 
throughout the wind and results in terminal wind velocities of the order of the escape speed at the base of 
the flow. 

The dynamics and radiation field of disk LDWs employed by Shlosman et al. (1985) and by Vitello & 
Shlosman (1988) were oversimplified. Both were approximated by a 1-D planar model allowing for divergence 
of the flow streamlines and geometrical dilution of the radiation field. Nineteen resonance lines in the range 
of 500— 1, 600 A were included in the calculation of the self-consistent radiation force. It was noted that disk 
LDWs are more restrictive than stellar winds and their development is strongly governed by the ionization 
structure in the wind. 

Subsequently, a variety of 2-D kinematical models for disk winds in CVs supplemented by a 3-D radiation 
transfer in the Sobolev approximation were explored (Shlosman & Vitello 1993; Vitello & Shlosman 1993). 
Calculations using an alternative Monte Carlo radiation transfer method, albeit with frozen-in ionization, 
gave similar results (Knigge, Woods, & Drew 1995). Constrained by synthetic line profiles and by calcu- 
lated effects of varying basic physical parameters, such as accretion and mass loss rates, temperature of the 
boundary layer, rotation, and inclination angle, the available phase space for wind solutions was sharply re- 
duced. Wind resonant scattering regions exhibiting a strongly bi-conical character regardless of the assumed 



velocity and radiation fields were identified and mapped. This allowed to match the observed line shapes 
from a number of CVs and to put forward a number of predictions, which were verified in high-resolution 
HST observations (Shlosman, Vitello, & Mauche 1996; Mauche et al. 1999). Most importantly, rotation was 
positively identified as the dominant factor shaping the UV line profiles in CVs confirming that the disk and 
not the white dwarf is the wind source. 

The above 1-D dynamical and 2-D kinematical modelings suffered from uniqueness problems which 
can be removed only by invoking the 2-D dynamics. Recent successful attempts by Proga, Stone, & Drew 
(1998; PSD hereafter) to model the 2-D time-dependent radiation hydrodynamics of disk LDWs was a major 
breakthrough in our understanding of this phenomenon. (The model of Murray & Chiang 1996 for CV winds, 
on the other hand, does not provide specifics for the wind geometry and hence we avoid discussing it here.) 
It confirmed basically that kinematical models of disk winds have sampled the correct parameter range and 
provided the scaling laws between different wind characteristics, e.g., between mass loss rate and accretion 
luminosity, and delineated the phase space for possible time-dependent solutions. A number of empirical 
relationships were put forward which require a physical explanation. 

In this paper we focus on the 2-D geometry of a disk LDW in the presence of a realistic radiation field 
in CVs. We analyze solutions of the wind Eulcr equation, emphasising differences in the solution topology 
with that of CAK stellar winds. In the subsequent Paper II we address issues related to the mass loss rates 
and velocity laws of CV winds. The possible contribution to wind-driving by magnetic stresses is ignored 
here (e.g., Blandford & Payne 1982; Pudritz & Norman 1986; Emmering, Blandford, & Shlosman 1992), as 
are jet-like outffows seen in other disk systems (Livio 1997). 

This paper is organized as follows. Section || reviews the relevant aspects of CAK theory for LDWs 
from O stars. Section y| addresses the 2-D geometry of disk LDWs, as well as the radiation field above the 
CV disk. Section H deals with an analytic solution for vertical winds above an isothermal disk, and Section 
analyzes the solution topology and flow geometry for tilted winds above a disk with a realistic temperature 
stratification. Section o compares our results with other models and observations, and Section n\ summarizes 
our work. 



2. CAK theory for O stars 

2.1. The stellar line force 

The CAK theory for line-driven winds from O stars forms the basis for our model of CV winds, and is 
therefore briefly summarized here. CAK assume a line-distribution function per unit v and k, from UV to 
IR, 

Ni.,.)^'-^(^y-\ (1) 

where v is the line frequency and k [cm^g~^] the mass absorption coefficient normalized to kq, where the 
latter refers roughly to the strongest driving line in the flow (Owocki et al. 1988). For the power exponent, 
< a < 1 holds, where the lower limit corresponds to purely optically thin lines, and the (unrealistic) upper 
limit to purely optically thick lines. Puis, Springmann, & Lennon (1999) derive a — 2/3 from Kramers' 
formula applied to resonance lines of hydrogenic ions. Similar values of a are obtained from detailed NLTE 
calculations for dense O supergiant winds (Pauldrach 1987; Pauldrach et al. 1994). On the other hand, 
for low-density winds, e.g., of B stars near-the-main sequence, a — 1/2 may be more appropriate (Puis, 
Springmann, & Owocki 1998). We shall, therefore, consider both cases a = 1/2 and 2/3 to study the effect 



gL = kq vth -; ~ \ c''^ A 7 ^7 "' (2) 



of a on the structure of disk winds. 

Using eq. (1), the CAK force from all lines can be written in a general way applicable for both geometries 
(e.g., Owocki & Puis 1996), 

r(a) 1 

\ — a (? 

by means of the Sobolev approximation (Sobolev 1957). V(a) is the complete Gamma function, dio is a solid 
angle centered about 7, and /^ is the frequency-integrated intensity in this direction. The line optical depth 
in direction 7 is given by 

"■'-7-V(7-v)' ^'^ 

with gas density p, and 7 • V(7 • v) being the gradient along 7 of the velocity component in direction 7. 
Note that Kg Wth is independent of the ion thermal speed Wth, and so is the line force. Assuming spherical 
symmetry, and adopting the 'radial streaming' approximation of CAK, i.e., t^ ^ Tr, eq. (0) simplifies to 

, ,!_„ r(a) F fdv/drY ,^, 

with frequency-integrated, radial flux F. 
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2.2. Stellar Euler equation 

For an isothermal, spherically symmetric stellar wind, the stationary Euler equation in dimensionless 
form can be written as, 

42\ 4A2 

1~^)W' = -1- — + FW-, (5) 

where, after CAK, we introduced a radial coordinate U = —B.^,/r, with i?* being the stellar radius. The 
sound speed. A, and the flow speed, V = VW, are normahzed to the photospheric escape speed from the 
reduced stellar mass, M(l — F), where F is the Eddington factor. The normalized wind acceleration is given 
by W = dW/dU = r'^vv' /GM{\ — F). Note that the gravitational acceleration is normalized to —1, whereas 
CAK normalize it to —1/2. The constant E in (^ is given by 
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where L is the stellar luminosity and M is the mass loss rate, M — inr^pv. Global solutions to (||) exist only 
above a certain, critical Ec, called an eigenvalue of the problem, i.e., below a maximum allowable mass loss 
rate. In the zero-sound speed limit, A = 0, the differential equation (|5|) separates into an algebraic equation, 

P = T^' -I- 1 - EW'" = 0, (7) 



and a trivial differential equation W' — const, which leads to the CAK velocity law, v == I'ocV'l ^ R*/f) 
where Woo is the flow terminal velocity. The Euler equation in the form (7) is particularly simple and its 
terms have a straightforward physical meaning, namely inertia, gravity and line force. 



2.3. The stellar wind topology: critical point of the flow 

We now consider solutions to eq. (ra) with finite A. According to CAK, for sufficiently large values of 
E there are two solutions in the supersonic regime, W > A^, termed 'shallow' (small W) and 'steep' (large 
W) solutions, whereas in the subsonic, photospheric regime, W < A"^, only the shallow solutions exist. On 
the other hand, only the steep solutions reach infinity. Namely, the term —AA^/U (the thermal pressure 
force due to geometrical expansion) becomes infinite for r ^ oo, and must be balanced by W — > oo along 
the branch of steep solutions. CAK concluded therefore that the true, unique wind solution has to switch 
from the shallow to the steep branch at a 'critical' point (see Fig. nl). Of course, a non-zero pressure term 
at infinity is unphysical because it requires an infinite amount of energy in the flow and is purely a result of 
the imposed isothermal conditions in the wind. 

The subsonic region has an extent of a few percent of the stellar radius for O star winds, while the 
pressure force —AA'^/U becomes important only beyond a few 100 i?*. In the intermediate regime, i.e., almost 
everywhere, the simplified eq. (H) with A = holds to a good approximation, and therefore W = const. 
This happens since both gravity and line force are ex r~^. Solution curves in the {W , U) plane are essentially 
straight lines, bending over due to thermal pressure both at C/ = — 1 and 0. As a result, the critical point is 
a saddle point in the {W , U) plane (Bjorkman 1995). The usual mathematical definition of a critical point 
refers, however, to the {W, U) plane or a topologically equivalent plane. 

For A = 0, the critical point lies on an extended ridge, and its position becomes ill-defined (Fig. [|a). 
In this limit, every point of the critical solution is a critical point. For A > 0, however, CAK find the 
critical point to lie at Tc = |-R* (Fig- |l|b). Inclusion of further correction terms to the line force, especially 
due to the finite size of the stellar disk, breaks the r~^ dependence of the line force, pushes the two almost 
degenerate critical solutions W'{U) apart, and shifts the critical point towards the sonic point (Fig. Hr). 
Pauldrach et al. (1986) ffiid then re ^ 1.1 i?*. 

A critical point is the information barrier for LDWs, and plays a role similar to the sonic point in 
thermal winds or nozzle flows (Abbott 1980). How can the pressure mismatch of a shallow solution at 
r > 300 i?* be communicated upstream to the critical point at 1.1 i?*? We speculate that it is not really the 
outer boundary mismatch which forces the flow through a critical point. Instead, the truly distinguishing 
property of the critical solution should be its correspondence to the maximum mass loss rate in the wind. 
Work is underway to identify the feedback mechanism between the wind and the photosphere which drives 
the wind from any shallow to the unique critical solution. This issue will be addressed elsewhere. In the 
present paper, we assum,e that the true disk wind solution is the one with the maximal allowable mass loss 
rate. Additional justification comes from the fact that only the shallow solutions for the disk wind connect to 
the disk photosphere. However, their terminal speeds are much smaller than the white dwarf escape speed, 
in sharp contrast to observed CV winds. The solution we are searching for should therefore switch to the 
steep branch (with large Voo) at some critical point, i.e., should be the solution of maximum mass loss rate. 

The flow critical point (subscript 'c') is defined by the singularity condition, dP/dW^ = (i.e., merging 
of a shallow and a steep solution). Together with the Euler equation, this implies 

W^ = -^, (8a) 

1 — a 

Ec = — -, — ^—r, — . (8b) 

a"(l-a)i-" ^ ' 

The eigenvalue Ec determines the maximum mass loss rate, and W^ determines the terminal speed. They 
are further discussed in Paper II. Furthermore, from the Euler equation, P = 0, also dP/dU = must 



hold everywhere. This leads to the regularity condition, dP/dUc + W^ dP/dWc = (if W^" < oo), which 
determines the position of the critical point. 



3. Disk wind geometry and radiation field in CVs 

3.1. Flow geometry, gravity and centrifugal force 

The central assumption throughout this paper is that the helical streamline of a fluid parcel in the 
wind is contained within a straight cone. While this is certainly an idealization, and a major restriction 
of this model, justifaction comes firstly from the related kinematical model of Shlosman & Vitello (1993); 
and secondly from the numerical 2-D hydrodynamic simulations of PSD, which showed that the escaping 
mass-loss carrying streamlines are well approximated by straight lines in the (r, z) plane, r and z being 
cylindrical coordinates. 

We denote the angle between the wind cone and the radial direction f in the disk plane by A. This angle 
is calculated self-consistently using the Euler equation, and is not assumed a priori. The footpoint radius of 
a streamline in the disk is Tq, and x is the distance along the cone (cf. Fig. 0). We search for the solutions 
of the Euler equation in the [A(ro), x] plane, for a streamline starting at arbitrary rg. The A(ro) dependence 
leads to the appearance of a new eigenvalue problem for the disk wind, and derivation of this function is the 
focus of the present paper. 

Since LDWs are highly supersonic, we neglect the pressure forces, and furthermore assume that the 
azimuthal velocity is determined by angular momentum conservation above the disk, and by Keplerian 
rotation in the disk plane. The tilt angle A has to be a monotonically decreasing function of tq to avoid 
streamline crossing, which would violate the assumption of a pressureless gas. The only remaining velocity 
component is v^, which points upwards along the cone. The dynamical problem has therefore been reduced 
to solving the Euler equation for Vx- 

In a frame rotating with the angular velocity J7 of a fluid parcel positioned at radius vector r, there are 
three fictitious forces (e.g., Binney & Tremaine 1987). The Coriolis acceleration, —2 $7 x r, has no component 
along X, and the same is true for the inertial force of rotation, —flxr. We introduce the effective gravity 
function, which is the component of gravity minus centrifugal force along the straight line in direction x, 
and is given by — (G Afwd/^o) .9(^i ^)' with M^d being the mass of the white dwarf, and 

(X \\^ -^ + '^Qs A cos A 

^^ ' ' (l + 2XcosA + X2)3/2 (l + XcosA)3' ^''> 

and X = x/tq. In the following, all lengths written in capital letters are normalized to the footpoint radius, 
To, in the disk. Close to the disk, X -^ 1, and g :^ X, while for X ^ 1, 5 ~ X~^. For a vertical ray, A = 90°, 
g has its maximum at X = l/\/2, while for a horizontal ray, A = 0, the maximum is at X = 1/2. 

Equation (|9|) shows an important difference between stellar and disk winds. The stellar gravity is always 
decreasing with distance, while for the disk an effective gravity 'hill' must be overcome before the wind can 
escape. This effect of disk LDWs will be discussed in Section |[ 



3.2. Radiation field above the disk 



Next, we evaluate the line force in eq. (0). Besides the initial growth of effective gravity with height, 
the opacity-weighted flux integral is the central property which distinguishes disk winds from stellar winds. 
Pereyra, Kallman, & Blondin (1997) give an analytical approximation for this integral above an isothermal 
disk. Unfortunately, an error was introduced with a change of integration variable, which led to an artificial, 
linear dependence of the vertical flux on z, even as 2; ^ 0. PSD solve this integral numerically (cf. Icke 
1980), using approximately 2,000 Gaussian integration points. 

In the general spirit of the radial streaming approximation of CAK, we replace the integral in eq. (y) 
by f ~"F, thus introducing an equivalent optical depth, f. We flrst calculate the frequency-integrated flux 
dF(r, z) at location (r, z) from a flat ring of radius q, radial width dq, and isotropic intensity I{q, 0), 



dF{r, z) 
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dF, 
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(10) 



where 

B = {r^ + z^ - q^f + Az\\ (11) 

For an isothermal disk with isotropic intensity, we integrate eq. diffl) over g, to obtain the disk flux 



^(^'^)"vi 



tt/ 1 f I[r2 + z2 + g2)] 
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(12) 



where r^ is the outer disk radius. For the nonmagnetic systems considered here, we identify the inner 
disk radius with r^jjj the white dwarf radius. We do not include contributions to the radiative flux from 
the white dwarf and the boundary layer. Generally, of course, accretion disks are not isothermal. We, 

— 1/2 

therefore, consider two complementary cases with T(ro) oc 7'q (termed 'Newtonian' disk in the following) 



-3/4 



1/4 



and T(ro) ex rg '^''* (1 — ^r„d/ro j (Shakura & Sunyaev 1973; SHS hereafter). Observations show that 
the brightness temperature stratification of CV disks is consistent with both these distributions (Home & 
Stiening 1985; Home & Cook 1985; Rutten et al. 1993). 



For the Newtonian disk, we find 



F(r, z) = 7r/(r)r^ 



/ 3r^ 



2r^/B 
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dz — r 
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(13) 



where 



C = 



(z2 _ r^)q^ + (z2 + ^2y ^ (^2 ^ ^2)^ 



(14) 



The surface flux above the SHS disk can only be integrated numerically, from eq. ([l(]|). Yet, this has 
the advantage that f can be introduced for each ring individually. More specifically, f is calculated along 
the fiux direction of a given ring at the position of the wind parcel. If the flux ( [l3|) is used instead, f is 
calculated along the disk flux direction. Typical differences in the resulting value for the tilt angle Ac (see 
below) are 5° to 10° for the two approaches. Corrections due to the t~" weighting in the azimuthal integral 
are even smaller. 



Figures ^ and ^ show isocontours for the z- and r-components of the flux (|3|) above the Newtonian 
disk. For sufficiently large tilt angles, the flux along the streamlines has a maximum larger than 7r/(ro, 0) at 
some X . This is due to the increasing visibility of the inner, hot disk regions. We denote this regime, where 
the flux has a maximum, the 'panoramic' regime, to be distinguished from the planar 'disk' regime, where 
Fz — ttI, and the 'far fleld' regime, where F oc X^^. 

We now introduce the normalized flux, F(r, z) — F(r, z)/7r/(ro, 0), along a streamline, which is indepen- 
dent of disk luminosity. To quantify the flux increase above the disk. Fig. g| shows the projected, normalized 
flux _F!j; = X • F as function of X , for different footpoint radii tq and for both types of non-isothermal disks. 
Note that the initial increase of F^ with X due to the wind's exposure to the central disk region is rather 
mild, a factor of a few only, because the central region has a small area. In Section u we discuss how the 
maximum in F^ controls the base extent of the wind above the disk, as well as the height of the wind critical 
point. 

In deriving the above fluxes, the intensity was assumed to be isotropic. Using instead the Eddington 
limb darkening law, Ig = |/o (l + | cos6'), with polar angle 9, the vertical flux in the planar disk regime 
above an isothermal disk becomes larger by a factor of 8/7, i.e., limb darkening should not significantly affect 
the wind properties. However, limb darkening can be more important in the UV spectral regime due to the 
Wien part of the spectrum, and the correction factors could become somewhat larger there (Diaz, Wade, & 
Hubeny 1996). 



4. Vertical wind above an isothermal disk 

As an analytically tractable case, we consider first a vertical (or cylindrical) wind with A ~ 90'', or x = z, 
above an infinite, isothermal disk with a flux F = tt/ z. We again adopt the 'radial streaming' approximation 
in (0), i.e., f = r^. Note that r^ has no contributions from either azimuthal velocity gradients, dv^/dr, or 
from geometrical expansion terms, oc w^, the latter describing photon escape along the tangent to the helical 
streamline. 

The density p which enters r is replaced by introducing the mass loss rate from one side of a disk annulus, 
dM. Since the mass which streams upwards between two cylinders is conserved (using again coordinate x), 

dAf(ro) = 2TrrodroVxiro,x) p{ro,x). (15) 

For simplicity, we apply the zero sound speed limit, ^ = 0, for the rest of this paper, and neglect the force 
due to electron scattering because of small F above the geometrically-thin disk. The Euler equation becomes 

= P{W',X) = W' + g-EW'°', (16) 

where g = X/{1 + X'^fl'^ for A = 90°, and W = 2YdYldX. Here djdX = tq djdx, and V is the flow speed 
normalized to the local escape speed at the footpoint tq on the disk. Normalizing instead to the escape speed 
from the white dwarf leads to unwanted, explicit appearances of rg in the Euler equation. The constant 
E{ro) for a streamline starting at tq on the disk is deflned as (compare with eq. |6|) 

E(r ) = ^^^^ ( ^°^*'' V~° 27rrgF^(ro,z = 0)/c^ 
1-a V27rGMwd/ (m dAf (ro)/dro)" 

Similarly to the stellar case, eq. (Qq) for the disk wind has solutions only above a critical value i?c, 
the eigenvalue of the problem, i.e., below a maximum allowable mass loss rate. Unlike the point star case. 



however, P in eq. (|6|) is a function of X even when A = 0. As a result, the degeneracy in the position of 
the critical point does not exist here, and one has a well-defined critical point, irrespective of A. 

There exists an additional difference between the stellar and disk LDWs. In Fig. np, the finite cone 
correction factor causes the critical point in the stellar wind to move upstream, and, for vanishing sound 
speed, the critical point and the sonic point both are found in the stellar photosphere. For the disk case, 
however, only the sonic point falls towards the photosphere, whereas the critical point stays at finite height. 
Namely, from the regularity condition dP/dX — (since P does not depend on W) , the critical point of the 
disk wind lies at the location of maximal gravity, at Xc = l/v2- 

This explains why Vitello & Shlosman (1988) find no critical point in the disk regime, X <C 1, for a 
vertical wind with a constant ionization. The variable wind ionization introduces additional gradients into 
the driving force, shifting the critical point towards the disk photosphere. For the solution discussed here, 
vertical ionization gradients are not mandatory. 

The conditions P = and dP/dW = lead to 



W^ = -^^.gc, (18a) 

1 — a 

Ec = — TT^-T-, ffc"", (18b) 

a"(l -«)!-" '^ ^ ^ 

where gc = 2/(3 v3)- This defines the wind solution of maximmn allowable mass loss rate. The effective 
gravity hill imposes a 'bottleneck' on the flow, i.e., the maximum of g{X) defines the minimum, constant 
eigenvalue Ec, or the maximum allowable Af , for the critical solution which extends from the disk photosphere 
to arbitrary large X. Larger values of E correspond to shallow solutions and hence to smaller mass loss 



rates. Smaller values of E correspond to stalling wind solutions. Note that Ec in eq. (18b) is independent 
of ro, in accordance with eq. (|lq). 



5. Tilted disk v^rinds 

With all pre-requisites at hand, we can now solve the general eigenvalue problem for a tilted wind above 
a non-isothermal disk. The density p in (^ is replaced by the conserved mass-loss rate between two wind 
cones, 

XroidX/droY 



dM{ro) = 2'Krodro (1 + X cos A) 



1- 



smXvxiro,x) p{ro,x). (19) 



sin A 

The term (1 + XcosA) describes the density drop due to the increasing radius of the cone, and [1 — 
Xro (dA/(iro)/sinA] describes the density drop due to the geometrical divergence of neighboring cones. 
The factor sin A stems from the quenching of the flow at small A. 



5.1. Disk Euler equation 

The geometrical expansion term V7 in the directional derivative 7 • V(7 • v) has contributions from the 
azimuthal curvature of helical streamlines and from the cone divergence dX/drQ. Close to the disk, where the 
mass loss rate of the wind is established, both contributions are small. For azimuthal curvature terms, this is 
shown in Appendix A. With regard to cone divergence, the argument is a posteriori, i.e., we find below that 
dX/dro is small. Two neighboring wind rays launched at, e.g., ro ~ 5rwd intersect at a normalized distance 



10 



Xi ~ —10 below the disk. Generally, X^ is larger by a factor of 10 than Xc, the distance between the disk 
plane and the critical point. By analogy with spherically symmetric stellar winds, where 7 • V(7 • f Vr) = 
/^^ dvr/dr + (1 — /i^) Vj./r, with /i = 7 • f , the geometrical expansion term oc v^ V(7 • x) for disk winds should 
be ex Vx/[ro{X ~ X^)]. Whereas the geometrical expansion term for O star winds, (1 — fi^)vr/r, is of the 
same order as the gradient term, /i^ dvr/dr, it is much smaller for disk winds. On the other hand, far from 
the disk, the expansion term may become important. However, we find from solving the Euler equation that 
it has only a marginal influence on the terminal wind velocity. Azimuthal terms for helical streamlines are 
unimportant far from the disk, where the wind is essentially radial. We, therefore, neglect all geometrical 
expansion terms in the following. Appendix A also shows that gradients in the azimuthal velocity can be 
neglected in the line force. Finally, we assume that the gradient of Vx points in the x-direction. This is a 
reasonable assumption since the velocity gradients develop roughly in the flux direction, as is also shown 
below. The normalized Euler equation for a conical disk wind, and for vanishing sound speed, is then 

O^P{W',X)^W' +g-EfW'°', (20) 

with auxiliary function /, 



firo,X) 



(1 + X cos A) f sin A - Xtq -— 
\ dro 



l+2a 



dF^i'+^". (21) 



I-wd 



Here, dF = |dF| (see eq. |T^), and ^ is the cosine of the angle between dF and the wind cone. Again, 
W — 2V dV/dX), where the velocity V is normahzed to the local escape speed; i?(ro) is defined in eq. (p^). 
Note that the flux integral in (p^) introduces a further dependence of / on tq. Furthermore, due to the 
weighting with /i^" in the integral, the disk flux vector and the wind cone do not generally point in the 
same direction. For disk winds as considered here, good alignment between radiative flux and wind flow is 
expected, however. In cases where such an alignment is not possible, e.g., for atmospheres irradiated from 
above, ablation winds at large tilt angle with the radiative flux were recently suggested (Gayley, Owocki, & 
Cranmer 1999). 



5.2. Wind tilt angle as an eigenvalue and solution topology 

The critical point conditions for a specific streamline are, from eq. ( |2C| ) 

W^c = 1 3c, (22a) 

1 — a 

Ec - — 7 r. — ^^, (22b) 

a"(l-a)i-" /e ' ^ ^ 

= (1-a)^-^. (22c) 

9c /c 



The tilted disk wind is essentially a 2-D phenomenon, and hence we expect two eigenvalues of the Euler 
equation, with respect to E and A. Finding the critical solution of maximum mass loss at a given footpoint 



ro implies minimizing E in eq. (22t) with respect to the position of the 'critical' point (Xc, Ac). We show now 
that (ATcjAc) is a saddle point of g^""//. We consider first the X coordinate, and recall from the analysis 
of the vertical disk wind that the maximum of g^^" has determined the eigenvalue Ec- From eq. ( P2b| ), the 
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relevant function now is g^^" / f ■ This means that the maximum of g^^°^ j j with respect to X for a fixed 
A serves as a bottleneck of the flow, i.e., the stringentest condition on the wind between the photosphere 
and infinity. It, therefore, defines the maximum allowable mass loss rate. Next, we analyze the mass loss 
rate along a streamline by varying its tilt angle A. To obtain the maximum mass loss rate, we look for the 
minimum of 5^~"// as a function of A. This particular Ac plays the role of a second eigenvalue of the Euler 
equation, besides E^- Note that due to the dependency of / on rp, the wind tilt will change with tq. The 
eigenvalue E^ is thus given by 

1 5^"" 

Ec — T- n niin max — - — . (23) 

This is the definition of a saddle point of g^^"' / f- Isocontours of this function are shown in Fig. ^ for the 
SHS disk. The existence of the saddle point in g^^" / f underlines the 2-D nature of disk LDWs. Because 
the saddle point opens in the X-direction, the wind escapes to large X. 

Furthermore, the critical solution of maximum mass loss passes also through a saddle point of the Euler 
function P in the {W , X) plane, in complete analogy with O star winds. The regularity condition, eq. ( p2q ), 
determines the loci Xc of these critical points, shown by the heavy lines in Fig. y. On the left branch of these 
curves, which passes also through the saddle point of g^~°'/f (if the latter exists), lie critical points of the 
saddle or X-type. Here, W'{X) can switch from a shallow (small W') to a steep (large W) solution. On the 
other hand, the right branch of the regularity curves, which passes through the minimum of g^^"//, consists 
of critical points of the focal type (Holzer 1977; Mihalas & Mihalas 1984). They correspond to solutions 
which do not extend from the disk photosphere to large radii, and are ignored in our discussion. 

Figure R shows a good overall alignment of the wind ray of maximum mass loss rate with the radiative 
disk flux vector, at least up to the critical point. This is (i) because the eigenvalue E^ depends linearly on /, 
but only with a small power of 1 — a on 5, and (ii) because only / has a maximum as function of A, whereas 
g falls off monotonically. 



5.3. Inner and outer disk winds 

Up to this point we ignored the possibility of multiple saddle points of 5^^"//. We now address this 
issue. As shown in Fig. g, for tq <^ 4rwd the function g^~°'/f has only one saddle at a large height, e.g., 
Xc — 4:A for a = 2/3. However, for ro ^ 4 7-^(1, a second saddle exists at smaller Xc, which lies on a different 
branch of the regularity curve. We name these two types of saddles the high and low saddle, according to 
their height Xc above the disk. The effective gravity 'hill' separates the two saddle points. 

^From Fig. ^, the low saddle corresponds to a larger mass loss rate than the high saddle. For tq ^ 4rwd, 
the solution of maximum mass loss is therefore determined by the low saddle. For smaller tq, however, only 
the high saddle exists, and determines the wind solution then. These two cases deflne the outer and inner 
disk wind, respectively. Clearly, the assumption of straight streamlines is a severe one for the inner wind 
with high- lying critical points. 

The tilt angle of the outer wind is around 60°, namely Ac = 65° at rg = 4rwd, and 55° at 20rwd- This 
is largely independent of a. For the inner wind, the tilt is larger, Ac = 80° for a = 2/3 and 70° for a = 1/2. 
Furthermore, the critical point, Xc, for the inner wind is much higher above the disk than the critical point 
for the outer wind. As mentioned above, these critical points fall on the opposite slopes of the effective 
gravity hill. For the outer disk wind, the position Xc of the critical point is moving closer to the wind sonic 
point with increasing rg . The reason for this is the larger gradient of the disk radiative flux in the x-direction 



12 



for larger tq. As a result, the line force can balance gravity at lower X. 

Figure ^ shows critical wind solutions W'{X) above the SHS disk, for different tq. The decelerating 
solution branches, W < 0, are discussed in Appendix B. The critical point topology of Fig. pi can be compared 
with that of the CAK stellar wind flow in Fig. H. (Note, that W has a slightly different definition for the 
stellar and disk wind cases.) From Fig. pi we can also derive a condition for the existence of a stationary, 
outer wind solution, further clarifying the role of the effective gravity hill. The plus signs at the critical 
points in the figure indicate where the Euler function P > 0, i.e., where drag forces (gravity and inertia) 
overcome the driving forces (line and centrifugal force), and correspondingly for the nfinus signs. Hence, 
d'^P/dXl < at the low saddle, or, using eq. (|2|), (1 - a) g'^/g^ < /"//c (resp. '>' at the high saddle). This 
means, the maximum of / must be sufficiently broad to allow for a stationary solution with a low saddle. 
Note that the critical point for a vertical wind above an infinite, isothermal disk, where f — 1, corresponds 
necessarily to a low saddle. 

In order to fully understand the geometry of disk LDWs, we consider also the transition region between 
the inner and outer winds. As discussed above, the low saddle does not exist below tq <^ 4rwd- Fig. 
shows g^^"/ f in the neighborhood of this footpoint radius. At tq — 4rwd, only the high saddle exists and 
determines the wind solution. At tq — 4.03 Twd, an inner regularity curve of elliptical shape has formed, but 
not yet the low saddle point oi g^^" / f . The mass loss rate is maximal at the smallest A along the curve, i.e., 
at its lower tip, which determines the wind solution in this transition regime. Finally, by rg = 4.15 r„d7 a 
low saddle has formed at Ac — 65°. The wind tilt stays at this value Ac, which corresponds to the maximum 
mass loss rate. In total, the wind tilt switches continuously from the high to the low saddle over a narrow 
range of 0.1 rwd in the footpoint radius. 



5.4. Overall disk wind geometry 

Table 1 lists important parameters of the wind above SHS and Newtonian disks, i.e., the tilt angle. A, 
the normalized mass loss rate from a disk annulus, m, and the critical point location, X^- The mass loss 
rate m is normalized to a vertical wind above an isothermal disk. The shallow maxima of the function F^ 
in Fig. ^ are responsible for m — 0{1). Implications of these mass loss rates are discussed in Paper II. From 
the table, one finds the ray dispersion in the outer wind, at intermediate footpoint radii from 4 to lOr^d, 

^--^ <^^' 

dro Twd 

Further in or out the ray dispersion is even smaller. Actually, since dX/dro enters the Euler equation (|20|), 
the full wind problem can be solved only iteratively. However, the dependence of the eigenvalues E^ and Ac 
on dX/dro is weak, and we assume throughout that eq. (p3) holds. 

The overall geometry of the disk wind is shown in Fig. |l^. For a = 2/3, the critical points are at 
Xc ~ 10 — 20rwd for the inner wind, then move towards the disk photosphere and stay at Xc — r^d, 
independent of footpoint radius tq in the outer wind. For a — 1/2, on the other hand, the critical points 
lie somewhat higher for the outer wind, at x^ — 2r^vd, but again independent of radius. While the division 
between the inner and outer wind persists (namely high-lying vs. low-lying saddle, or critical points on 
opposing sides of the gravity hill), the transition in A between the two regions is smooth for a = 1/2, and 
the inner tilt reaches a maximum of A = 70°. 

The innermost disk region, 1 — 2 rwd, is left out in Fig. nu. The details of the disk wind and its very 
existence here are subject to great uncertainties in the radiation field which depends on the properties of the 
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transition layer and the white dwarf itself. The outer boundary of the disk LDW is set by the radius where 
the disk temperature falls below 10 K and the wind driving becomes inefficient, in analogy with stellar 
winds (Abbott 1982; Kudritzki et al. 1998). For the SHS disk with L^ = 10 L©, this should happen around 
30 7-wd- 



6. Discussion 

We compare here our theoretical model of LDWs from accretion disks in CVs with those available in the 
literature, kinematical and dynamical ones. We ignore the radial wind models, with the white dwarf being 
the wind base, because they are in a clear contradiction with current observations (e.g., review by Mauche 
& Raymond 1997). An alternative source of gas is the disk itself. Kinematical models which account for 
this source of material subject to the line-driving force successfully explained the observed bipolarity of the 
outflow, and reproduced the inclination-dependent line profiles (Shlosman & Vitello 1993). Their weak point 
was the absence of a unique solution. The 1-D dynamical models in a simplified disk radiation field revealed 
some major differences between the stellar and disk winds, e.g., the bi-polarity and the existence of a gravity 
hill (Vitello & Shlosman 1988). 

More sophisticated 2-D kinematical models, supplemented with a 3-D radiation transfer in Sobolev 
approximation, showed the importance of rotation in shaping the lines (Vitello & Shlosman 1993; Shlos- 
man et al. 1996). Finally, the 2-D hydrodynamical model of a disk wind in a realistic radiation field and 
with the line-force parameterized by the CAK approximation has addressed the issue of flow streamlines and 
mass loss rates in the wind (PSD). Our comparison, therefore, is focused on these models. 

Vitello & Shlosman (1993) set up a kinematical disk wind model assuming straight flow lines in order 
to fit the C IV P Cygni line profiles of three CVs observed with the lUE. The fit parameters included the 
inner and outer terminating radius of the wind base, and the corresponding tilt angles of the wind cone. The 
best fit appeared to be indifferent to the mass loss rate, within the range of 10^^ to 10^^ of the accretion 
rate. In the present work, which accounts for wind dynamics, we find lower mass loss rates more justified 
and discuss various implications of these rates on the wind models in Paper II. The tilt of the innermost 
wind cone in Vitello & Shlosman was rather steep, A — 80°, while at the outer disk edge A — 25°. A similar 
work by Knigge et al. (1995), but using Monte Carlo radiation transfer in the wind, gave similar results. 
In the present work, the tilt angle A is calculated self-consistently from the Euler equation, resulting in a 
similar inner tilt as found from kinematical models, while the outer tilt differs by a factor of 2 between the 
two approaches. 

The most advanced numerical modeling of CV winds from the SHS disk so far was performed by PSD, 
using the time-dependent Zeuss 2-D code. We find a number of similarities between their and our results, 
but differences exist as well. Our comparison with PSD is limited to their models 2-5, i.e., without a central 
luminous star. These models agree with ours on the overall wind geometry. This includes the streamline 
shapes and the run of the wind opening angle with height. The streamlines in PSD appear to form straight 
lines in the (r, z) plane, in striking similarity with the previous kinematical models. In addition, the change 
in the wind opening angle with distance from the rotation axis seems to be weak in PSD. The mass loss rates 
are consistent between both models, and so arc the wind optical depths, which can approach unity even for 
very strong resonance lines (Paper II). 

While PSD also find two markedly distinct flow regions, the inner and outer, their inner wind, at 
To <j 4r„dj appears as the sole and only outflow. The outer disk region, at radii ^ 4rwd, exhibits a time- 
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dependent irregular flow, resulting in no mass loss. In contrast, in our model, mass loss from the SHS 
disk is dominated by the inner wind and the innermost part of the outer wind, as is discussed in Paper II. 
Interestingly, our outer wind seems to be more robust than the inner wind. For the latter, the balance of 
driving and drag forces which leads to a high saddle on the far side of the gravity hill is rather a delicate 
one. Setting, for example, the centrifugal force arbitrarily to zero causes the high saddle solution to vanish, 
whereas the low saddle remains almost unchanged. 

PSD suggest that the irregular behavior of the outer flow is a consequence of the different X-dependence 
of gravity and disk flux, with the gravity preventing the wind from developing. We find instead that at radii 
ro ^ 4r„d, where a low saddle exists, the fast increase in the projected disk flux, Fx{X), results in a 
sufficiently strong growth of the line force which drives the wind past the gravity hill. For the inner wind 
regions, on the other hand, where no lower saddle branch exists, the wind indeed must overcome the gravity 
barrier without the appropriate radiation flux increase with X. 

Based on our analysis of the disk wind we may provide some insight into the erratic behavior of the 
flow PSD observe at larger radii. It is possible that mass overloading of the fiow causes spatial and temporal 
fluctuations of the streamline divergence dX/dvo, wherefore the wind stalls on characteristic length scales 
of the vertical gravity, tq. We conjecture that fiuctuations on such scales may be self-amplifying, or, in 
other words, result in a locally converging flow with dX/drQ > 0. However, if such an instability exists, it 
is expected to be most pronounced in the inner disk regions. Namely, with increasing tq, the flux increases 
faster with X, which moves the critical point closer to the sonic point. Mass overloading seems therefore 
less likely in the outer wind regions, in contrast to the findings of PSD. A linear stability analysis of our 
stationary wind model with respect to harmonic perturbations of dX/dro would also be interesting, but is 
beyond the scope of the present work. 

Furthermore, we cannot confirm the dependence of A on the disk luminosity as in the PSD model. 
We find that the eigenvalue Ac for each streamline is determined from the positions of the saddle points 
of the function g^~°' /f- Both g and / are independent of the disk luminosity, / specifically so because it 
is normalized to the flux at the streamline footpoint (eq. |2|). Therefore, Ac depends only on the radial 
temperature stratification in the disk. 

One important issue neglected in our modeling is the saturation of the line force at some value when 
all the driving lines become optically thin. If this thick-to-thin transition occurs before the flow reaches its 
critical point, the wind solution is lost, since the drag forces overcome the driving forces. However, this 
still leaves the possibility that a more complicated wind dynamics is established, where the decelerating flow 
at some larger radius starts again to accelerate (i.e., jumps from a WL to a W[_ solution). We leave this 
question open for future scrutiny, and note here that the mass loss rates derived from the present eigenvalues 
E are upper limits. 

The present work is based on the CAK theory for stellar winds. Over the years, questions have been 
raised concerning the physical meaning of the CAK critical point (Thomas 1973; Lucy 1975; Cannon & 
Thomas 1977; Abbott 1980; Owocki & Rybicki 1986; Poe, Owocki, & Castor 1990). Most interesting 
for the present context is the inclusion of higher order corrections to the diffuse line force in the Sobolev 
approximation, which shift the critical point still closer to the sonic point (Owocki & Puis 1998; see also 
Fig. ^). This proximity of the sonic and critical points may not be coincidental, and one can speculate 
whether or not the sonic point determines the mass loss rate instead of the critical point. However, we find 
for the disk wind model that the sonic and critical points occasionally lie far apart, e.g., for a vertical wind 
above an isothermal disk, or a tilted wind close to the rotation axis ('inner wind'). 
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These fundamental issues impair our understanding of LDWs from stars and disks, and therefore must 
be addressed in the future. 



7. Summary 

We discuss an analytical model for 2-D stationary winds from accretion disks in cataclysmic variable 
stars. All parameters chosen are typical for high-accretion rate disks in novalike CVs. We solve the Euler 
equation for the wind, accounting for a realistic radiation field above the disk, which drives the wind by 
means of radiation pressure in spectral lines. Some key assumptions are that each helical streamline lies on 
a straight cone; that the driving line force can be parameterized according to CAK theory; and that the 
thermal gas pressure in the supersonic wind can be neglected. Our results can be summarized as follows. 

The disk wind solutions are characterized by two eigenvalues, the mass loss rate and the flow tilt angle, 
Ac, with the disk. The additional eigenvalue Ac for each streamline reflects the 2-D nature of the model. We 
find that the wind exhibits a clear bi-conical geometry with a small ray dispersion. Specifically, two regions, 
inner and outer, can be distinguished in the wind, launched from within and outside 4rwd, respectively. The 
tilt angle for the outer wind is A ^ 60° with the disk. At these angles, the wind flow and radiative disk flux 
are well aligned. For the inner wind, the tilt angle is larger, up to 80°. We emphasize that the disk wind 
tilt angle (i.e., the wind coUimation) depends solely upon the radial temperature stratiflcation in the disk, 
unless there is an additional degree of freedom such as central luminosity associated with nuclear burning 
on the surface of the white dwarf. 

A major distinction between stellar and disk winds is the existence of maxima in both the gravity and 
the disk flux along each streamline. The latter flux maximum appears to be a crucial factor in allowing the 
wind to pass over the gravity 'hill'. The flux increase is more pronounced further away from the rotation 
axis. As a result, the critical point of the outer wind lies close to the disk photosphere and to the sonic point. 
In fact, it lies upstream of the top of the gravity hill, and this proximity of the critical and sonic points is 
typical of LDWs from O stars as well. On the other hand, for the inner wind, the increase in radiation flux 
with height is smaller, and the critical point lies far away from the sonic point, beyond the gravity hill. 

Comparing our analytical models with the 2-D numerical simulations of Proga et al. (1998), we find 
an overall good agreement in the streamline shape, tilt angle, and mass loss rate, but our wind baseline is 
wider. 

We are grateful to Jon Bjorkman, Rolf Kudritzki, Chris Mauche, Norman Murray, Stan Owocki, Joachim 
Puis, and Peter Vitello for numerous discussions on various aspects of line-driven winds. I. S. acknowledges 
hospitality of the IGPP/LLNL and its Director Charles Alcock, where this work was initiated. This work 
was supported in part by NASA grants NAG5-3841 and WKU-522762-98-06, and HST AR-07982.01-96A. 



A. Line force due to gradients in the azimuthal velocity 

We estimate here the importance of azimuthal velocity terms for the line force in x-direction. Assuming 
Keplerian rotation within the disk, and angular momentum conservation above the disk, one has 



dv^/dz 1 1 1 vW 

dv^/dx ~ ""71 tan A 1 -F X cos A ~W 
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dv^/dr _ 1 1-XcosA Vw 

T2 "FTTT; (Al) 



dv^/dx V2 {l + Xcos\f W" 

__Hh_ _ \/2 %/W 

Here, the singularity of tan^^ A at A = is a result of neglecting the pressure terms in the Euler equation. 



Note that dv^/dr changes sign at X = 1/cosA. From (Al), gradients in v^ are comparable to gradients 
in Vx when ^/W /W ^ 1. The main question is for their influence on the mass loss rate. Because, in the 
CAK model, M is determined by the conditions at the critical point, we calculate \/W /W at the latter. We 
consider first a vertical wind from an isothermal disk. Since W grows monotonically up to and somewhat 
beyond the critical point (see Fig. |), and because W = ^ W'dX, one has y/Wc/W^ < (3/2)3/4^1 _ a/a - 1. 
Here, Xc — l/\/2 and eq. ( |l8a|) were used. Alternatively, the critical points of the outer wind above a non- 



isothermal disk typically lie close to the disk, hence g{X) ~ X. Using eq. (22a), y/X^JWl ~ y/l — a/a ~ 1. 
Both disk cases give, therefore, essentially the same result. We conclude that w^-terms can be important 
everywhere between the disk photosphere and the critical point, and hence can modify M. 

To find their effect on M, we include u^-terms in the evaluation of the line force, eq. (g), in an ap- 
proximate manner. Only the disk regime is considered, in which case the radiation intensity is roughly 
isotropic and the radiation flux has a z-component only. The azimuthal part of the solid angle integral in 
(0) is approximated by a 4-point quadrature at angles mr/2 with f, where n = 0,1,2,3. This leads to a 
correction factor of the approximate form |;(2 -I- |1 -f 5*1" + |1 — 5*1") to the line force, EfW". Here, S 
is a linear combination of the expressions in eq. (|Al|), with coefficients < 1 from angle integration. In the 



disk regime, S'c '^ 1 from eq. (Al). This coincides with the borderline between an increase and a decrease 
in M due to inclusion of u^-terms, which lies at Sc = 1.25 for a — 1/2 and at Sc — 1.18 for a = 2/3. A 
detailed, numerical calculation of the above angle integral is required to decide, which of both cases actually 
occurs. Since, however, Sc is close to unity, the influence on the mass loss rate is limited to a 10% effect. 
We, therefore, neglect w^-terms in calculating the line force. 



B. Disk wind deceleration 

In Figure O, isocontours which cross through the low saddle point loop into one another at some larger 
height, Xd- At X > Xd, one has E > Ec from the flgure, i.e., the allowed maximum mass loss rate in this 
region is smaller than that at the saddle. At these distances, inertia and gravity overcome the line force plus 
centrifugal force, and the wind decelerates, W < 0. As is shown in Paper II, the wind speed always exceeds 
the local escape speed at X^, which implies that the decelerating wind reaches inflnity at a positive speed. 

Due to the deceleration, the velocity law becomes non-monotonic, and the line transfer is no longer 
purely local, because global couplings occur between distant resonance locations. We neglect these couplings 
here, and simply replace W'"' in the line force by [VF' |". For a wind ray launched at ro = 5 r^d, Fig. || shows 
that a single, decelerating branch, WL < 0, accompanies the critical, accelerating solution W[_ of maximum 
mass loss rate. It is suggestive that at X^ the solution curve jumps from the W'^ to the W'_ branch, and 
extends thereupon to infinity. 

The discontinuity in W' introduces a kink in the velocity law. Such kinks propagate at sound speed 
(Courant & Friedrichs 1948; actually, for LDWs, at some modified, radiative-acoustic speed, see Abbott 
1980 and Cranmer & Owocki 1996) and are therefore inconsistent with the assumption of stationarity. It 
seems plausible, however, that the discontinuity in W' is an artifact of the Sobolev approximation, since 
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the latter becomes invalid at small dv/dx, i.e., as W — > 0. An exact line transfer should instead give 
a smooth transition from W!^_ to W'_- We find inc 
dW^ _/dX(i —f — oo, e.g., in the top panel of Fig. ||. 



a smooth transition from W!^_ to W'_- We find indeed cases of 'almost' smooth transitions, where both 
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Table 1. Tilt angle Ac, normalized mass loss rate TOc, and position X^ of the critical 
point for different disk wind models. Underlined numbers indicate the transition from the 

inner to the outer wind. 

SHS Disk Newtonian Disk 

a = 2/3 a == 1/2 a = 2/3 a = 1/2 

Tn/r^d Ac[°] rhc X^ K rfic Xc Ac rfic X^ Xc "^c Xc 
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Fig. 1. — Saddle point topology of stellar CAK winds in the {W',U) plane, for (a) a point star and zero 
sound speed; (b) a point star and finite sound speed; (c) an extended star and finite sound speed. Filled 
dots mark the sonic points (at U ~ —1) and the critical points. 
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Fig. 2. — Adopted flow geometry for a CV disk wind. The streamlines are helical lines, and are assumed to 
lie on straight cones. 
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Fig. 3. — Isocontours of the frequency-integrated, vertical flux component Fz above a Newtonian disk. The 
disk extends from 1 to 30rwd- Normahzation is /(rg = Sr^d, z = 0) = 1. 
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Fig. 4. — Same as Fig. 0, including normalization, here for the radial flux component F,-. Dotted lines 
indicate an inward flux. 
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Fig. 5. — Normalized, projected flux Fx{roTX) for SHS (top panel) and Newtonian (bottom panel) disks, 
at footpoint radii ro — 20, 10, 5, and 2rwd (top to bottom curves). The tilt angle with the disk plane is 90° 
(full lines) and 60° (dashed hues). 
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Fig. 6. — Isocontours oi g^~°'/f, normalized to (2/3-v/3)^~", over the {X,X) plane. Footpoints of the wind 
are at tq = 3, 5, and ISrwd (left to right panel). The isocontours have logarithmic spacing. The temperature 
stratification is that of the SHS disk {r^ — SOr^jd, Oi — 2/3). Heavy hues are solutions to the regularity 
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Fig. 7. — Radiative flux vector above the SHS disk (ta = 30rwd)- The straight line indicates a wind cone 
at eigenvalue A — 61° (cf. Table 1). The '+' sign marks the critical point. 
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Fig. 8. — Critical wind solutions W'{X') at zero sound speed, above the SHS disk [r^ — 30rwd, a = 2/3). 
At ro — B^wd: high saddle; at rg = 5 and 15rwd: ^ow saddle. The branches W < are decelerating wind 
solutions, and are discussed in Appendix B. '+' and '— ' signs refer to the sign of the Euler function P in 
eq. (|20|). For comparison with the critical point topology of a spherically symmetric, stellar wind, see Fig. |lj. 
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Fig. 9.— Isocontours of g^""//, normalized to (2/3\/3)i""; from 0.855 to 0.885 in steps of 2 x lO"^ ^t 
the low saddle; from 1.05 to 1.15 in steps of 10~^ at the high saddle. At tq — 4, 4.03, and 4.15 rwd (left to 
right panel). Heavy lines are solutions to the regularity condition. Circles mark the critical wind solution 
of maximum mass loss rate. Within the footpoint range from tq = 4 to 4.15 rwd, the wind switches from 
Ac — 80° to Ac — 65° via a growing regularity curve of ellipsoidal shape. 
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Fig. 10. — Wind geometry above the SHS disk {r^ — SOr^d) according to Table 1, for a ~ 2/3 (left panel) 
and 1/2 (right panel). Black regions indicate the accelerating LDW, and thin white lines show a few wind 
cones. Grey areas indicate decelerating wind (Appendix B), for a ray dispersion dX/dro = — l°/rwd- 'I' 
and 'O' mark the inner and outer wind, 'T' is the transition region. Heavy, white lines are locations of flow 
critical points. The innermost region of the disk, at tq < 2r„d7 is not treated due to uncertainties in the 
radiation field. 



